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also the effects of quantum corrections. 
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Nowadays there is considerable renewed interest in the transition from 
classical mechanics to quantum mechanics, a powerful motivation behind that 
being the problem of the so-called quantum chaos (see Ozorio de Almeida 
1990, Gutzwiller 1990, Casati and Chirikov 1995). An important aspect is 
represented by the semiclassical quantization formula of the (regular) energy 
levels for quasi-integrable systems (Maslov and Fredoriuk 1981, Rau 1992, 
Braun 1993), the so-called torus quantization, initiated by Einstein (1917) 
and completed by Maslov (1972). 

It has been recently shown (Gram and Paul 1987, Degli Esposti, Gram 
and Herczynski 1991) that, for perturbed non-resonant harmonic oscilla- 
tors, the algorithm of classical perturbation theory can be used to formulate 
the quantum mechanical perturbation theory as the semiclassically quan- 
tized classical perturbation theory equipped with the quantum corrections 
in powers of H "correcting" the classical Hamiltonian that appears in the 
classical algorithm. In effect, one can explicitly calculate corrections to the 
Einstein-Brillouin-Keller (EBK) quantization of the classical tori (Maslov 
and Fredouriuk 1981). For example, the quantum corrections of the one- 
dimensional x A perturbed harmonic oscillator have been studied in great 
details by Alvarez, Gram and Silverstone (1988). 

Examples of rather detailed studies of semiclassical approximations, their 
resummations and of the Birkkoff-Gustavson normal forms can be found in 
Ali, Wood and Devitt (1986) and Ali and Devitt (1989). Another example of 
a rather complete semiclassical analysis of a one-dimensional system, namely 
the quartic oscillator, has been published by Voros (1983). 

The aim of this paper is to extend previous studies to a two-dimensional 
system, which is more interesting because it is nonintegrable and thus generic. 
The integrable systems are rather exceptional in the sense that they are 
typically isolated points in the functional space of Hamiltonians and their 
measure is zero in this space. If we randomly choose a system realized in 
nature, the probability is one that the system is nonintegrable (see Robnik 
1995). 

The model is given by two non-resonant oscillators coupled by a nonlinear 
quartic interaction of strength g (Pullen and Edmonds 1981): 

H =^1 + ^ + ^1 + ^1) +9^1 (1) 
Note that a similar Hamiltonian has been obtained for the Yang-Mills-Higgs 
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quantum mechanics (Savvidy 1984, Salasnich 1995). 

Through the canonical transformation in action-angle variables (see Born 
1927 and Dittrich and Reuter 1992): 

q k = \f2hcos6 k , p k = ^J21k~sm6 k , A; =1,2 (2) 
the Hamiltonian can be written: 

H = H (I 1 ,I 2 )+gV(I 1 ,I 2 ,e 1 ,e 2 ), (3) 

where: 

H (h,I 2 ) = wJi + W2I2, (4) 

V(h, I 2 , 6 U 6 2 ) = AIJ 2 cos 2 0x cos 2 2 . (5) 

Following the classical perturbation theory (see Born 1927, Dittrich and 
Reuter 1992), we search for a canonical transformation 1 2 , 9 1: 2 ) — > 
(Ii, I 2 , 0i, 9 2 ) to obtain a new Hamiltonian that depends only on the new 
action variables up to the second order in a power series of g: 

H(h,I 2 ) = H Q {h,I 2 ) + gH 1 (h, I 2 )+g 2 H 2 (h,I 2 ). (6) 

The generator S of the canonical transformation is supposed to be capable 
of expansion as a power series in g of the form: 

S(h, I 2 , 0x, e 2 ) = hO, + J 2 2 + gS 1 (I 1 , 1 2 , 0x, 2 ) + g 2 S 2 (h, J 2 , 0!, 2 ), (7) 

and to satisfy the equations: 

dS ~ dS x 2 dS 2 

Ik = W k =Ik + 9 W k +9 W k > (8) 

~ dS n dS 1 2 dS 2 
oh dl k dl k 

From the Hamilton-Jacobi equation: 

(10) 
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we have a number of differential equations that result on equating the coef- 
ficients of the powers of g: 

H {h, k) = H {h, ~h) = wi/i + cu 2 I 2 , (11) 

Htih, h) = (u,^ + ^2^) + ^(/l, J 2| 01, 2 ), (12) 

fwf fx / 03, dS* dVdS x dVdS, 

The unknown functions Hi, Si, H 2 and S 2 may be determined by averaging 
over the time variation of the unperturbed motion. 
At the first order in g we obtain: 

Hi(h,h) = -L r r d6id6 2 V(IiJ 2 ,6i,6 2 ) = hh, (14) 
4n z Jo Jo 



and 



Si{Ii, I 2 , 6i, 9 2 ) = --Jih[— sin 29i + — sin 29 2 + 
+ sin 2{6i - 9 2 ) + sin 2(0! + 6 2 )}. (15) 

UJ X - UJ 2 UJi + ^ 2 

At the second order in g we have: 

fr/r n 1 /" 27r /' 27r j/i j/i ,9V dSx dVdSi, 

^4U[4(A + A ) _(^ + M ] . (16) 

The approximate Hamiltonian (6) depends only on the actions f], so that 
by an application of the EBK rule: 

I k = (n k + ^)h, (17) 



3 Note that the integrable approximate Hamiltonian (6) could be obtained alternatively 
as the Birkhoff-Gustavson normal form, which is a purely algebraic method of calculating 
the action variables order by order for perturbed harmonic oscillators with polynomial 
perturbations (Robnik 1984, Kaluza and Robnik 1992, Kaluza 1993, Robnik 1993). 
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to equations (11), (14) and (16), we obtain a semiclassical analytical formula 
of the energy levels. 

Now we show how to connect this semiclassical formula with the usual 
quantum perturbation theory. In quantum mechanics the generalized coor- 
dinates satisfy the usual commutation rules [qk,Pi] = ih~8kh with k, I = 1,2. 
Introducing the creation and destruction operators: 

1 l 

Ofc = -/=(?* + m), al = — =(q k - ip k ), (18) 



<2h V2h 
the quantum Hamiltonian can be written: 

H = H + gV, (19) 

where: 

1 1 

H — huj^d^ai + -) + hu 2 (a% a 2 + -), (20) 

h 2 

V = -{a l + a+)\a 2 + a+) 2 . (21) 

At this point we mention the general problem of quantization (Robnik 
1984, Robnik 1988, Abraham and Marsden 1978): There is no unique quan- 
tization prescription. If one chooses certain quantization procedure, this 
will not commute in general with the classical canonical transformations. 
However, the quantization does commute with the linear classical canonical 
transformations. This is exactly the approach implemented in our present 
case (18), and thus our quantization is equivalent to the coordinate space 
quantization (the latter one yields the right quantum mechanics, whose re- 
sults agree with experiments). 

If \ri\n 2 > is the basis of the occupation numbers of the two harmonic 
oscillators, the matrix elements are: 

< n'^lHoln^ >= fibril + ^) + ^("l + j^l^n/nana, (22) 



and: 



n 



2 



< n x n 2 \V\n x n 2 >= — [y^iK - l)5 n ; ni _ 2 +V(ni + 1)K + 2)5 n ' ini+2 +(2rai+l)£ n / n Jx 



x [\lMn2 ~ l)5 <n2 _ 2 + J(n 2 + l)(n 2 + 2)5 n , n2+2 + (2n 2 + l)^ n J. (23) 
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The Rayleigh-Schrodinger perturbation theory (see Messiah 1962) up to the 
second order gives us: 

E(nih, n 2 h) = E (nih, n 2 Ti) + gE^n-Ji, n 2 Ti) + g 2 E 2 (n 1 h, n 2 h), (24) 

where: ^ 

E (n 1 h,n 2 k) = 7j[wi(ni + -) + u 2 (n 1 + -)], (25) 



E 1 (nih,n 2 h) =< riin 2 \V\nin 2 >, (26) 

| < n 1 n 2 \V\nin 2 > | 2 
h[ui(ni - rai) + uo 2 (n 2 - n 2 ) 



, -.s x ^ I < n\n'r,\V\nin 2 > I 2 



"1"2 



(n 1 ,n 2 )^(ni,n 2 ) 

We obtain immediately: 



E^m^n^) = a 2 (m + ^)(n 2 + ^), (28) 

and after some calculations: 

ft 3 m(ni-l)n 2 (n2-l) K + 1)K + 2)(n 2 + l)(ra 2 + 2) 

E 2 {nxTi,n 2 Ti) = — ■ ■ h 

32 a;i + u; 2 ^ + tf 2 

| wi(wi - l)(w 2 + l)(n 2 + 2) _ (m + + 2)n 2 (n 2 - 1) | 

a;i — uj 2 LO\ — uj 2 

| mjm - l)(2n 2 + l) 2 _ (m + l)(ni + 2)(2n 2 + l) 2 | 

(2wi + l) 2 n 2 (n 2 - 1) _ (2m + l) 2 (n 2 + l)(n 2 + 2) 

CJ 2 CJ 2 

The zero and first order quantum terms coincide with the semi-classical ones: 
E (mh, n 2 K) = H (( ni + ^)H, (n 2 + ^)h), (30) 

E 1 (n 1 h,r^h) = H 1 {(n 1 + ^)h,(r^ + ^)h), (31) 
and the second order quantum term can be written: 

E 2 ( ni h, n 2 h) = H 2 {{m + ^)h, (n 2 + ^)h) + h 2 Q 2 {{m + ^)h, (n 2 + ^)h), (32) 



where: 

2( / 1 ,/ 2) ._l [ M + M ] . (33) 
32 wi - u 2 lo\ + uj 2 

The quantum series rearranges directly into the classical canonical perturba- 
tion series plus quantum corrections proportional to successive powers of h 
(see Degli Esposti, Graffi and Herczynski 1991). 

The term h 2 Q 2 represents the quantum corrections to the EBK quantiza- 
tion up to the second order of perturbation theory. These quantum correc- 
tions depend linearly on quantum numbers. It is in contrast with first order 
results in a previous paper of Robnik (1984) in which the spectra differ only 
by an additive constant independent of the quantum numbers (actions). 

We compute the "exact" levels with a numerical diagonalization of the 
truncated matrix of the Hamiltonian (19) in the basis of the unperturbed 
oscillators (see Graffi, Manfredi and Salasnich 1994). The numerical energy 
levels depend on the dimension of the truncated matrix: We compute the 
numerical levels in double precision increasing the matrix dimension until 
the first 100 levels converge within 8 digits (matrix dimension 1225 x 1225). 

Then we compare the "exact" levels to the levels of the quantum per- 
turbation theory and to the semiclassical results. A very good agreement is 
observed for the lowest energy levels (see Table 1). 

In table 2 we show the error in units of the mean level spacing D between 
the "exact" levels and the levels obtained with semiclassical and quantum 
perturbation theory. 

We observe that the algorithm provided by the appropriate semiclassical 
quantization is comparable to the algorithm provided by ordinary quantum 
perturbation theory but the quantal corrections do not always increase the 
accuracy. Thus for some of the calculated levels the semiclassical quanti- 
zation gives better results than the quantum perturbation theory. In fact, 
Rayleigh-Schrodinger perturbation theory diverges for any value of g, but it 
has been proved to be Borel summable to the exact energy levels only for 
some special systems like the anharmonic oscillators with / degrees of free- 
dom with a polynomial perturbing potential which is asymptotically positive 
definite (Graffi, Grecchi and Simon 1970, Reed and Simon 1978, Caliceti, 
Graffi and Maioli 1980). 

In Figure 1 we plot the error, in units of the mean level spacing D, 
between the "exact" levels and the semiclassical levels. The figure shows 
some systematics but the prediction of individual levels by EBK quantization 
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worsens as the quantum number increases, contrary to the naive expectation 
(see also Prosen and Robnik 1993). As is seen, by decreasing H from 1 to 
1CT 1 the quality of the approximation improves considerably but only for 
sufficiently low levels (see also Table 3). 

If h, no matter how small, is kept fixed, the EBK quantization (torus 
quantization) of the individual levels is only a first order approximation of 
an expansion in h. Therefore, the accuracy of the approximation decreases 
for higher levels. To get a good agreement it is necessary, as is well known, to 
implement the classical limit, i.e. h — > and ni,n,2 — > oo, while at the same 
time keeping the actions I\ = {n\ + 1/2)% and I2 = (712 + l/2)h constant 
(Graffi, Manfredi and Salasnich 1994). 

In Figure 2 we plot the error, in units of the mean level spacing D, between 
semiclassical levels and levels obtained with quantum perturbation theory. 
Also in this case the error increases for higher levels. 

In conclusion, we have examined the transition between the classical 
and the quantum mechanics for a two-dimensional, nonintegrable and non- 
resonant system. Up to the second order of perturbation theory we have 
decomposed the quantum description into the classical description (i.e. the 
leading semiclassical term) plus quantum corrections which depend linearly 
on the quantum numbers. The analytical energy levels are in good agreement 
with the "exact" numerical ones, the semiclassical quantization is compara- 
ble to the quantum perturbation theory, and for some levels the semiclassical 
quantization gives better results than quantum perturbation theory. 

Finally we note that the extraction of quantum corrections for resonant 
systems is a more intricate procedure; some initial results for perturbed res- 
onant oscillators can be found in Graffi (1992). 

The present paper supports the following general conclusions. 

(cl) The quantum corrections, defined as the difference between the semi- 
classical and the quantum perturbation terms of the same order, as derived 
by Degli Esposti, Graffi and Herczynski (1991), certainly vanish in the limit 
h -> 0. 

(c2) The above conclusion, however, changes drastically when the quan- 
tum corrections are measured in natural units of the mean level spacing D, 
which scales as D ~ fi/, where / is the number of freedoms: The quantum 
corrections and their low fr-power terms can even diverge or remain constant. 
This is implicit in Degli Esposti, Graffi and Herczynski (1991). Probably the 
same conclusion applies when exact levels are compared to the semiclassi- 
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cal ones, supported by the results in Prosen and Robnik (1993) and by the 
results of the present paper. 

(c3) Since both the classical and quantum perturbation series typically 
diverge and thus do not necessarily describe the exact levels, not even after 
a certain resummation (except for some important notable exceptions previ- 
ously described (Gram, Grecchi and Simon 1970)), it is important to compare 
the semiclassical approximation (and the quantal perturbation results) with 
the exact spectra, which in general is impossible, since we generally do not 
have explicit solutions of the Schrodinger problem in a closed form. There- 
fore, we stress the importance of specific case studies like the present one, in 
order to get a better understanding of the quality of semiclassical approxi- 
mations. 

Acknowledgements 
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Table 1: Comparison between "exact" levels and levels obtained by pertur- 
bation theories. First 25 levels. E ex are "exact" levels, E sc are semiclassical 
levels, and E qp are levels obtained with quantum perturbation theory, h — 1, 
g = 10 -1 , LO\ — 1 and w 2 = \/2- 
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E ex_ E scy D \ E ex_ E qpy D 



1.0611359-10 3 
1.5578579-KT 2 
8.1338054-10- 3 
4.5591835-10" 2 
3.6215890-10- 2 
9.3476856-10- 3 
9.9898852-10" 2 

0.1227176 
7.0703819-10- 2 

0.1829406 
2.9902905-10~ 2 
0.2880960 
0.2973495 
0.2030921 
0.1350395 
4.085474M0- 2 
0.5432989 
0.4450069 
0.4778005 
0.2858459 
0.8948155 
0.3193286 
0.5559518 
0.1297049 
0.9053394 



1.1284242-10 
7.1859419-10 
9.0260478-10 
2.6613854-10 
4.2791300-10 
2.2781115-10 
7.0337765-10 
0.1187100 
9.2249520-10 
0.1428019 
5.0089385-10 
0.2735053 
0.2466222 
0.2140520 
0.1715501 
0.1029161 
0.5181223 
0.3836938 
0.4781800 
0.3117708 
0.8590558 
0.2980870 
0.5567808 
0.2067311 
0.8951331 



-2 



Table 2: The error measured in units of the mean level spacing D for the 
first 25 levels. D is calculated for the lowest 100 levels. E ex are "exact" 
levels, E sc are semiclassical levels, and E qp are levels obtained with quantum 
perturbation theory, k — 1, g — 10 _1 , ui\ — 1 and uj 2 = a/2- 
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Table 3: The error measured in units of the mean level spacing D between 
"exact" levels and semiclassical levels. First 25 levels. D is calculated for 
the lowest 100 levels, g = 10" 1 , — 1 and lu 2 = \/2. 
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Figure Captions 



Figure 1: The error measured in units of the mean level spacing between 
"exact" levels and semiclassical levels. First 100 levels, g = 10 _1 , — 1 
and uj 2 = (a) h = 1, (b) h = 10" 1 . 

Figure 2: The error measured in units of the mean level spacing between 
semiclassical levels and levels obtained with quantum perturbation theory. 
First 100 levels, g = 10" 1 , u x = 1 and uo 2 = (a) h = 1, (b) h = 10 _1 . 
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